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ABSTRACT 

The origin of radio emission from plerions is considered. Recent observations suggest 
that radio emitting electrons are presently accelerated rather than having been in- 
jected at early stages of the plerion evolution. The observed flat spectra without a 
low frequency cutoff imply an acceleration mechanism that raises the average particle 
energy by orders of magnitude but leaves most of the particles at the energy less than 
about few hundred MeV. It is suggested that annihilation of the alternating mag- 
netic field at the pulsar wind termination shock provides the necessary mechanism. 
Toroidal stripes of opposite magnetic polarity are formed in the wind emanated from 
an obliquely rotating pulsar magnetosphere (the striped wind). At the termination 
shock, the flow compresses and the magnetic field annihilates by driven reconnection. 
Jump conditions are obtained for the shock in a striped wind. It is shown that post- 
shock MHD parameters of the flow are the same as if the energy of alternating field 
has already been converted into the plasma energy upstream the shock. Therefore 
the available estimates of the ratio of the Poynting flux to the matter energy flux, cr, 
should be attributed not to the total upstream Poynting flux but only to that associ- 
ated with the average magnetic field. A simple model for the particle acceleration in 
the shocked striped wind is presented. 

Key words: acceleration of particles - magnetic fields - MHD - shock waves - 
pulsars: general - supernova remnants 



1 INTRODUCTION 

Most of the pulsar spin-down power is carried away by a 
relativistic, magnetized wind. The pulsar wind injects this 
energy into the surrounding nebula in the form of relativis- 
tic electron-positron pairs and magnetic fields therefore such 
nebulae, or plerions, emit synchrotron radiation from the ra- 
dio to the gamma-ray band. The most famous and well stud- 
ied example of the plerion is the Crab Nebula; the spectrum 
of this source is measured from about 10 MHz to dozens 
TeV. According to MHD models (Rees & Gunn 1974; Ken- 
nel & Coroniti 1984; Emmering & Chevalier 1987; Begelman 
& Li 1992), the pulsar wind terminates in a standing shock 
at a radius defined by the condition that the confining pres- 
sure balances the momentum flux of the wind. In the Crab 
case, the shock radius was estimated to be 3 x 10 17 cm in 
excellent agreement with observations; according to Chan- 
dra results (Weisskopf et al. 2000), the radius of the shock in 
the equatorial plane is 4 x 10 17 cm. The observed brightness 
and the spectral index distributions are generally consistent 
with the assumption that the relativistic particles are accel- 
erated at the termination shock and then fill in the nebula, 
spending the acquired energy on synchrotron emission and 
pdV work. 

The generic observational feature of plerions is a flat ra- 



dio spectrum; the spectral flux may be presented as a power 
law function of the frequency, T v oc v~ a , with the spec- 
tral index a — — 0.3. At high frequencies the spectrum 
steepens and the typical spectral slope in the X-ray band is 
q <; 1. The overall spectrum of the Crab may be described 
as a broken power law with the spectral breaks around 10 
Hz, fewx 10 15 Hz and around 100 keV. The synchrotron life- 
time of the radio emitting electrons (and positrons, below by 
electrons I mean both electrons and positrons) significantly 
exceeds the plerion age therefore one can not exclude a pri- 
ori that they were injected at the very early stage of the 
plerion evolution (Kennel & Coroniti 1984; Atoyan 1999). 
In this case the overall spectrum depends on history of the 
nebula. However the spectral break at 10 13 Hz may be sim- 
ply accounted for the synchrotron burn off effect assuming 
that particles emitting from the radio to the optical bands 
are injected more or less homogeneously in time with the 
single power law energy distribution. This view is strongly 
supported by Gallant & Tuffs (1999, 2002) who found that 
the infra-red spectral index in the central parts of the Crab 
is close to that in the radio, and gradually steepens outward. 
Recent observations of wisps in the radio band (Bietenholz 
& Kronberg 1992; Bietenholtz, Frail & Hester 2001) suggest 
unambiguously that the radio emitting electrons are accel- 
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erated now in the same region as the ones responsible for 
the optical to X-ray emission. 

If the radio emitting electrons have been injecting into 
the Crab Nebula till the present time, the injection rate of 
electrons should be about 10 40 — 10 41 s _1 . It is interest- 
ing that the observed pulsed optical emission from the Crab 
pulsar suggests that about the same amount of electrons is 
ejected from the pulsar magnetosphere (Shklovsky 1970) so 
the pulsar does supply the necessary amount of particles. 
The observed spectral slope in the radio band, a = 0.26, 
implies the energy distribution of the injected electrons of 
the form N(E) oc E~ 1,5 . In this case, most of the particles 
find themselves at the low energy end of the distribution 
whereas particles at the upper end of the distribution dom- 
inate the energy density of the plasma. Taking into account 
that no sign of a low frequency cut-off is observed in the 
Crab spectrum down to about 10 MHz whereas the high 
frequency break lies in the ultra-violet band (recall that the 
break at about 10 13 Hz is attributed to synchrotron cooling 
but not to the injected energy distribution), one concludes 
that the above distribution extends from E min < 100 MeV 
to E max ~ 10 6 MeV. At E > E max the distribution be- 
comes steeper; the spectral slope in the X-ray band, a = 1.1, 
corresponds, with account for the synchrotron burn off ef- 
fect, to N(E) oc E~ 2 ' 2 . The distribution further steepens at 
E ~ 10 9 MeV as it follows from the gamma-ray spectrum 
of the Nebula. Thus the injected electrons have a very wide 
energy distribution, their number density being dominated 
by low energy electrons whereas the plasma energy density 
being dominated by TeV electrons. 

The above considerations place severe limits on the pul- 
sar wind parameters and possible mechanisms of the par- 
ticle acceleration at the termination shock. According to 
the widespread view, the pulsars emit Poynting-dominated 
winds however the electro-magnetic energy is efficiently 
transferred to the plasma flow such that the magnetization 
parameter a, defined as the ratio of the Poynting flux to 
the kinetic energy flux, is already very small when the flow 
enters the termination shock. The mechanisms of such an en- 
ergy transfer still remain unclear (the so called cr-problem) 
however dynamics of the flow in the Nebula suggests that 
the magnetic pressure is small just beyond the termina- 
tion shock, which implies that the Poynting flux just up- 
stream the shock is very small (Rees & Gunn 1974; Kundt 
& Krotscheck 1980; Kennel & Coroniti 1984; Emmering & 
Chevalier 1987; Begelman & Li 1992). In this case it is ki- 
netic energy of the upstream flow that converts into the 
energy of accelerated particles when the plasma flow is ran- 
domized at the shock. Then the characteristic downstream 
"temperature" is about the upstream particle kinetic en- 
ergy, T ~ mcT m , so the average particle energy does not 
vary considerably across the shock. A high energy tail may 
be formed in the particle energy distribution (the parti- 
cle acceleration at the relativistic shocks is considered by 
Hoshino et al. (1992); Gallant & Arons (1994); Bednarz & 
Ostrowski (1998); Gallant & Achterberg (1999); Kirk et al. 
(2000); Achterberg et al. (2001)) however this tail merges, 
at its low energy end, with the quasi-thermal distribution at 
E ~ T ~ mc r„ . Therefore if the pulsar spin-down power, 
L sd , is converted into the kinetic energy of the wind, the 
available upper limit on the low-frequency break in the Crab 



spectrum suggests that the wind Lorentz-factor, F w , does 
not exceed few hundred. 

On the other hand, the observed flat spectrum may 
be formed only if the energy per electron in the wind is 
much larger than m e c 2 F w . Gallant et al.(2002), modifying 
the original idea of Hoshino et al. (1992) and Gallant & 
Arons (1994), suggested that the wind is loaded by ions; in 
this case the radio emitting electrons are accelerated by res- 
onant absorption of ion cyclotron waves collectively emit- 
ted at the shock front. The necessary ion injection rate, 
~ L sd /(m p c 2 T w ) ~ 10 39 s _1 , vastly exceeds the fiducial 
Goldreich- Julian elementary charge loss rate, ~ 3x 10 34 s _1 . 
Although one can not exclude by observations that pulsars 
emit the required amount of ions, the available pulsar mod- 
els do not assume an ion outflow with the rate exceeding 
the Goldreich- Julian charge loss rate (Cheng & Ruderman 
1980; Arons 1983). 

Here I consider an alternative possibility, which does not 
imply a radical modification of the basic pulsar model. In the 
equatorial belt of the wind from an obliquely rotating pul- 
sar magnetosphere, the sign of the magnetic field alternates 
with the pulsar period forming stripes of opposite magnetic 
polarity (Michel 1971, 1982; Coroniti 1990; Bogovalov 1999). 
Observations of X-ray tori around pulsars (Weisskopf et al. 
2000; Helfand, Gotthelf & Halpern 2001; Gaensler, Pivo- 
varoff & Garmire 2001; Pavlov et al. 2001; Gaensler et al. 
2002; Lu et al. 2002) suggest that it is in the equatorial belt 
where most of the wind energy is transported; theoretical 
models (e.g., Bogovalov's (1999) solution for the oblique split 
monopole magnetosphere) support this conclusion. There- 
fore the fate of the striped wind is of special interest. In the 
striped wind, the Poynting flux converts into the particle 
energy flux when the oppositely directed magnetic fields an- 
nihilate (Coroniti 1990; Lyubarsky & Kirk 2001, henceforth 
LK; Lyutikov 2002; Kirk & Skjaeraasen 2003). Until now this 
dissipation mechanism was considered only in the unshocked 
wind. It was found that the flow acceleration in the course 
of the energy release dilates the dissipation timescale so that 
the wind may enter the termination shock before the alter- 
nating field annihilates completely. In this case driven anni- 
hilation of the magnetic field at the shock may provide the 
energy necessary to form the flat particle distribution. On 
the other hand, the formed distribution may extend down 
to low enough energy because the kinetic energy of the flow 
in the striped wind is lower than the total energy. The aim 
of this research is to consider properties of the termination 
shock in the striped wind. 

It will be shown that the alternating field completely 
annihilates at the shock so that the downstream parameters 
of the flow are the same as if the field has already anni- 
hilated upstream the shock. Therefore the available limits 
on the upstream magnetization parameter, a, should be at- 
tributed not to the total Poynting flux but to the Poynt- 
ing flux associated with the averaged magnetic field (see 
also Rees & Gunn 1974; Kundt & Krotscheck 1980). The 
upstream flow may be Poynting dominated provided most 
of the Poynting flux is transferred by alternating magnetic 
field. On the other hand, driven reconnection of the magnetic 
field within the shock radically alter the particle accelera- 
tion process. It follows from both analytical and numerical 
studies that particles are readily accelerated in the course of 
the magnetic field reconnection to form the energy distribu- 
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tion with the power-law index j3 ~ 1 (Romanova & Lovelace 
1992; Zenitani & Hoshino 2001; Larrabee, Lovelace & Ro- 
manova 2002). Electrons with such a distribution emit syn- 
chrotron radiation with a flat spectrum therefore Romanova 
& Lovelace (1992) and Birk, Crusius-Watzel & Lesch (2001) 
suggested that reconnection plays a crucial role in flat spec- 
trum extragalactic radio sources. One can naturally assume 
that radio emission of plerions is generated by electrons ac- 
celerated in the course of reconnection of the alternating 
magnetic field at the pulsar wind termination shock. Then 
the steeper high-energy spectrum may be attributed to the 
Fermi acceleration of particles preaccelerated in the recon- 
nection process. A simple model for the particle accelera- 
tion in the shocked striped wind is presented here. It will 
be shown that the minimal energy of the power law en- 
ergy distribution may be low in this case, even less then 
the upstream kinetic energy, whereas the energy density of 
the plasma will be dominated by high-energy electrons in 
agreement with the observed spectra of plerions. 

The article is organized as follows. Jump conditions for 
a shock in the striped wind are obtained in Sect. 2. Making 
use of these jump conditions I demonstrate in sect. 3 that 
the alternating magnetic field dissipates completely at the 
pulsar wind termination shock. In sect. 4 I discuss particle 
acceleration by driven reconnection. Particle acceleration at 
the shock in the striped wind is analyzed in sect. 5. The 
obtained results are summarized in sect. 6. 



2 JUMP CONDITIONS FOR THE SHOCK IN 
A STRIPED WIND 

Let us find jump conditions for a shock in a flow with an al- 
ternating magnetic field. The plasma upstream the shock is 
assumed to be cold everywhere with the exception of narrow 
current sheets separating stripes with opposite magnetic po- 
larity. The pressure balance implies that the magnetic field 
strength in adjacent stripes differs only by sign but not by 
the absolute value; on the other hand the width of stripes 
with opposite polarity may not be the same so that the 
average magnetic field may be nonzero (Fig. 1). One gener- 
ally finds jump conditions from the conservation laws and 
some prescription for the magnetic flux passing the shock. 
In the standard MHD shock, the magnetic flux is conserved; 
in the case under consideration one should choose more gen- 
eral prescription to take into account possible annihilation 
of the magnetic field at the shock. 

Outside the shock, the magnetic field is frozen into the 
plasma, E + (l/c)v x B = 0. In the case of interest the 
magnetic field lies in the plane of the shock. In the frame 
where the flow velocity is perpendicular to the shock plane 
E = (v/c)B. Then the Faraday law may be written as 

dt c dx 
Together with the continuity equation 



this implies 



d d 

_r„ + -r„„ = o, 



I! 






Figure 1. The magnetic field in the striped wind. The average 
field is shown by dashed line. 



where V and n are the Lorentz factor and the proper plasma 
density, correspondingly, and b is a constant for each fluid 
element. The energy and momentum fluxes are 

a n 2 . EB 
b = wl v A — - — c; 

47T 



wY 2 {v/c) 



■P + 



B 



871 



where p and w are the gas pressure and specific enthalpy, 
correspondingly. Taking into account the above considera- 
tions, one can write 



S = WvT 2 



WT 2 (v/c) 2 +V, 



where the effective pressure and enthalpy are 



W = w + 



(2) 



(3) 



8tt ' 4tt ' 

so the flow may be described by the hydrodynamical equa- 
tions. 

In the standard MHD theory, the magnetic flux con- 
serves and the flow is considered as homogeneous both up- 
stream and downstream the shock so that b is a global con- 
stant. In our case 6 is an alternating function of the la- 
grangian coordinate however only b 2 enters the conservation 
laws therefore the only essential difference from the stan- 
dard MHD theory is that the modulus b may decrease in 
the shock because of the field annihilation. Let us introduce 
a parameter 



6? 



(4) 



(1) 



to measure decreasing in the magnetic flux across the shock. 
The indexes 1 and 2 are referred to quantities upstream and 
downstream the shock, correspondingly. 

One can express the downstream parameters via the 
upstream ones making use of the conservation of the particle, 
energy and momentum fluxes across the shock. Radiation 
losses may be safely neglected at the shock width scale. In 
the striped wind, one should take the fluxes averaged over 
the wave period. Let us assume, for the sake of simplicity, 
that the current sheets between the stripes are so narrow 
that one can neglect the contribution of the plasma within 
the sheets into the conserving fluxes. Outside the sheets, 
the magnetic field and plasma density are constant across 
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a stripe and the adjacent stripes are differ only by the sign 
of the magnetic field and by the width. Taking into account 
that the energy and momentum fluxes are proportional to 
the magnetic field squared, one can write the conservation 
laws like in the case of homogeneous magnetic field: 



Wiwir? = w 2 « 2 rl; 



mviiv^c) 2 +Vi = mrl(v 2 /cy + r 2 . 



(5) 
(6) 



(7) 



The only formal difference from the case of the standard 
MHD shock is in the additional parameter r) < 1; in the 
standard shock ij = 1. 

One can resolve the above equations as follows. The 
plasma downstream the shock is relativistically hot, W2 = 
4p2, therefore one can easily eliminate w and p from Eq.(3) 
and express Vi via W2 and ri2. Substituting the obtained 
relation into Eq.(7), one then eliminates W2 making use of 
Eq.(6) and eliminates 712 making use of Eq.(5). Taking into 
account that the plasma upstream the shock is cold, p\ = 
and Mi = nm e c 2 , one gets the resulting equation for the 
downstream velocity in the form 



Vl — V2 



r?«i 



■q v 2 F 2 v 1 T 1 



where 



4 v\Yl Auv 2 V 2 ' 



4:nm e c 2 



(8) 



(9) 



is the magnetization parameter. This equation is exact. In 
case r? 3> Max(u, 1) one can take vi = c. Then Eq.(8) 
reduces to 

(3V2-C) = ^(C + V2). 

\ a ) V2 

Now the downstream parameters may be found explicitly: 

v 2 = I (i + x + \/i + i 4 x + x 2 ) ; (10) 



n 2 = niTi 



^2[17-8 X -X 2 - (i + x)yr+14 X + X 2 ] 



i + x+ v /l + 14 x + : 



(11) 



^ = £ T^V /l7 - 8 ^-^ 2 -( 1 + ^v /l + 14 x + x 2 



6x 



1 + X + a/1 + 14X + X 2 



where 



X = 



l + o-' 



(12) 



(13) 



At 77 = 1 one recovers the downstream parameters for 
the relativistic MHD shock in the homogeneous medium 
(Kennel & Coroniti 1984; Appl & Camenzind 1988). In this 
case the quantity x is the ratio of the upstream Poynting 



flux to the total energy flux. In a striped wind, \ is deter- 
mined by the Poynting flux corresponding to the magnetic 
flux passed the shock. It will be shown in the next section 
that the alternating magnetic field annihilates completely at 
the pulsar wind termination shock; then x is simply the ratio 
of the Poynting flux associated with the average magnetic 
field, {Bi) 2 v/4tv, to the total energy flux. At the equator 
of the flow the average field is zero, a — 0, therefore one 
gets V2 = c/3 like at the nonmagnetized relativistic shock. 
When the average magnetic flux is large, a 1, r) — > 1, the 
downstream flow is relativistic with the Lorentz factor 



r 2 = 



1 + (l-Tj)ff' 



(14) 



3 THE TERMINATION SHOCK IN A 
STRIPED PULSAR WIND 

In a striped wind, the magnetic field forms toroidal stripes of 
opposite polarity, separated by current sheets (Michel 1971, 
1982; Bogovalov 1999). Such a structure arises in the equa- 
torial belt of the pulsar wind and may be considered as an 
entropy wave propagating through the wind. Usov (1975) 
and Michel (1982) noticed that the waves must decay at 
large distances, since the current required to sustain them 
falls off as r _1 whereas the density of available charge carri- 
ers in the wind decreases as r~ 2 . Hence the alternating mag- 
netic fields should eventually annihilate. It was shown in LK 
that the distance beyond which the available charge carriers 
are unable to maintain the necessary current exceeds the 
radius of the standing shock where the wind terminates so 
that only some fraction of the magnetic energy is converted 
into the particle energy before the plasma reaches this shock 
front. This fraction depends on the magnetic field reconnec- 
tion rate. A lower limit may be obtained assuming that the 
dissipation keeps the width of the current sheet to be equal 
to the particle Larmor radius, which is roughly the same 
condition as that the current velocity is equal to the speed 
of light (Corontiti 1990; Michel 1994; LK). Then about 10% 
of the Poynting flux dissipates before the wind reaches the 
termination shock. At higher reconnection rate the fraction 
of the dissipated energy is larger; assuming that the recon- 
nection velocity is as high as the sound velocity one can find 
that the alternating field annihilates completely before the 
flow enters the termination shock (Kirk & Skjasraasen 2003). 
Here I adopt the slow reconnection model by LK to estimate 
the flow parameters just upstream the termination shock. 

Neglecting the dependence of the wind parameters on 
latitude, one can express them only via the pulsar spin-down 
power, L a d, and the total amount of the ejected electrons, N. 
LK used instead the ratio of the gyro-frequency at the light 
cylinder to the angular velocity of the neutron star and the 
so called multiplicity coefficient arising in the models of the 
pair production in pulsar magnetospheres. These parameters 
may be expressed via L sd and N taking into account that 
th e magnetic field at the light cylinder is estimated as Bl = 
y L s d/(cRi), where R t — cP/2-k is the light cylinder radius, 
P the pulsar period. 

The characteristic distance beyond which the available 
charge carriers are unable to maintain the necessary current 
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(Eq.(14) in LK) may be written as 



sheets is 



Rma 



e 2 L s 



■Ri 



(15) 



1.9 x 



where e is the electron charge. For the Crab Rmax 
10 19 cm. Note that the dimensionless parameter L used by 
Kirk & Skjasraasen (2003) is proportional to (R ma x/Ri) 2 - 
In the Poynting dominated flows, dissipation of even a small 
fraction of the Poynting flux implies significant acceleration 
of the flow. While the dissipated energy is small, the Lorentz 
factor of the wind (Eq.(30) in LK) may be estimated as 



where 



r -ir 

Z V ■tirna. 



r - 

J- max — 



c 2 N 



(16) 



(17) 



is the Lorentz factor attained by the wind if all the spin- 
down power is converted into the kinetic energy of the 
plasma. Note that Y w /Y max is the fraction of the spin-down 
power transferred to the plasma therefore the magnetization 
parameter of the wind may be written as 



- 1. 



(18) 



The current sheet width is conveniently measured by a frac- 
tion A of a wavelength 2nRi occupied by the two current 
sheets; Eq.(31) in LK may be written as 



R 



Rn 



(19) 



In the Crab, the termination shock is observed at the dis- 
tance 4 x 10 17 cm (Weisskopf et al. 2000) therefore just up- 
stream the shock a = 13, A = 0.024; F w = 4.4 x 10 3 /iv4o, 
where A> 40 = iV/(10 40 s" 1 ). 

At the termination shock, the flow sharply decelerates 
so the plasma is compressed. The proper density grows 
Ti/r2 times. According to Eqs.(10) and (13), the less r) (i.e. 
the larger fraction of the magnetic flux dissipates at the 
shock) the less the downstream velocity so one can place a 
lower limit on the compression factor assuming r\ = 1. In 
this case r 2 = ^Ja; substituting the estimated above param- 
eters of the Crab wind, one gets Fi/r2 = V w j^fa ~ 1000. 
Such a huge compression factor implies significant heating 
of the plasma, especially within the current sheets. Let us 
show that after such a compression the particle Larmor ra- 
dius should exceed the wavelength so that the alternating 
magnetic field dissipates completely within the shock. 

The downstream temperature, which determines the 
Larmor radius, depends on the fraction of the alternating 
magnetic field annihilated at the shock. Let us first as- 
sume that this fraction remains small so that the struc- 
ture sketched in Fig. 1 preserves downstream the shock. The 
width of the current sheet in the wind frame cannot be less 
than the particle Larmor radius, 



S = T/eB', 



(20) 



where T is the plasma temperature within the sheet and 
B' the magnetic field in the wind frame of reference. The 
wavelength in the wind frame is 2nRiT therefore the mini- 
mum fraction of the wavelength occupied by the two current 



A = T/(nRiTeB') 



The temperature obeys the condition of the hydrostatic 
equilibrium of the plasma in the sheet, n^T = B' 2 /8tt, where 
rih is the number density of the plasma in the sheet. The 
magnetic field is frozen into the plasma outside the sheet 
therefore B' oc n c , where n c is the plasma number density 
outside the sheet. Now one can write (see Eq.(17) in LK) 

n c r 

The left-hand side of this relation is the fraction of parti- 
cles carried in the sheet; of course this fraction should be 
less than unity. In the unshocked wind rih/n c = 3 (LK); 
with the above estimate for the upstream A, one obtains 
that upstream the shock this fraction is about 0.1. It was 
demonstrated above that T decreases about 1000 times when 
the flow passes the shock, therefore downstream the shock 
Arih/ric ~ 100, which is impossible. This means that con- 
trary to the initial assumption, the fraction of the magnetic 
field dissipated within the shock is not small and r\ defined 
by Eq.(4) is not close to unity. 

In this case x (Eq.(13)) is also not close to unity and 
then, according to Eqs.(10) and (12), the downstream flow is 
non- or mildly relativistic, r 2 ~ 1, whereas the downstream 
temperature is about the particle kinetic energy upstream 
the shock, T 2 ~ Y w am e c 2 . With the above estimates for 
the upstream Lorentz factor and magnetization parameter, 
one can easily see that the Larmor radius downstream the 
shock, Vl = T/eB ~ 10 1 /(B-4N40) cm, vastly exceeds the 
wavelength, A = 2irRi — 10 cm, therefore the alternating 
magnetic field should completely annihilate at the shock. 

This means that the parameter \ should be determined 
as the ratio of the Poynting flux associated with the average 
magnetic field, c{B) 2 /47r, to the total energy flux. Therefore 
the downstream parameters are independent of whether the 
alternating magnetic field dissipated upstream the shock or 
just at the shock. The available upper limits on the magne- 
tization parameter in the pulsar wind (Kennel & Coroniti 
1984; Emmering & Chevalier 1987; Begelman & Li 1992) was 
found from the standard MHD shock conditions applied to 
the downstream parameters estimated from the analysis of 
the plasma dynamics in the nebula downstream the shock. 
In the case of the striped wind these upper limits should be 
attributed not to the total Poynting flux but to the Poynting 
flux associated with the average magnetic field (see also Rees 
& Gunn 1974; Kundt & Krotscheck 1980). The upstream 
flow may be Poynting dominated provided most of the elec- 
tromagnetic energy is carried by the alternating magnetic 
field, which annihilates at the termination shock. 



4 PARTICLE ACCELERATION BY DRIVEN 
RECONNECTION 

The particle acceleration by driven reconnection may be 
qualitatively described as follows. Let us consider a box with 
two stripes of the oppositely directed magnetic field and a 
current sheet between them (Fig. 2). When the box is com- 
pressed, the electric field E = — (l/c)v x B arises, which has 
the same sign in the domains of opposite magnetic polarity. 
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Figure 2. A compressing box with oppositely directed magnetic 
fields. Distribution of the magnetic field is shown at the bottom. 
The current sheet is shaded. The magnetic field is directed per- 
pendicularly to the plane of the figure. Directions of the velocities 
and the electric fields are shown by thick arrows. An unbounded 
particle trajectory near the field reversal is shown by dashed ar- 



Close enough to the zero line of the magnetic field, particles 
may move freely along this line and gain energy from the 
electric field. Of course the real picture should be much more 
complicated than the presented one-dimensional sketch. The 
magnetic reconnection might proceed in separate X-points 
but in average the process remains one-dimensional because 
the released energy is confined to a layer around the field 
reversal. 

Acceleration of relativistic electrons close to an X-point 
was considered by Romanova & Lovelace (1992); Zenitani & 
Hoshino (2001), Larrabee et al.(2002) who found a power- 
law particle distribution with the slope /3 ~ 1. Reconnection 
in a long current sheet at a time-scale large enough that 
particles pass many X-points has not been considered yet. 
Let us assume that a power-law distribution is formed in this 
case also and estimate the maximal energy particles attain 
when the magnetic field annihilates completely. 

Let the particle energy distribution within the current 
sheet be N(-y) = K^'* 3 at 1 < 7 < 7™ with /? ~ 1 - 
2. At such a distribution function, the particle density is 



dominated by low energy electrons, 

K 

rih 



(21) 



whereas the energy density is dominated by high-energy 
electrons, 



m e c 2 K 2-/3 



■P 



■In 



(22) 



Let us assume that the power law index (3 remains fixed in 
the course of compression and only K and 7 m vary. One 
can find these variations considering particles and energy 
balance within the box. 

The particle balance is written as 



n h 8 + n c (l — S) = n h0 6 Q + n c0 (lo — So), 



(23) 



where 5 is the sheet width, / the box width, nu and n c the 
particle number densities in the sheet and outside it, corre- 
spondingly, and the index is referred to the initial state. 
The energy balance implies that the variation of the total 
energy within the box is equal to the work done on the box 
by the outer pressure (= B 2 /8ir): 



d 



e5+(l- 5) 



Stt 



= - — dl 

8tt 



where e is the plasma energy density in the sheet. Taking 
into account that the magnetic field is frozen into the cold 
plasma, 



B = bn c 



(24) 



and that the plasma pressure in the sheet, p = e/3, is coun- 
terbalanced by the magnetic pressure, 

(25) 



e _ B_ 

3 ~ 8tt' 

one can write the energy balance equation in the form 



,, dn c dS 



(26) 



Let us assume that the sheet width is equal to the max- 
imal Larmor radius 

m e c 2 7 m 

The maximal Lorentz factor, j m , is found from Eqs. (21), 
(22), (24) and (25) as 



87r(/3 — 1) meC 2 rih 



Now one can write 



S = 5i 



\n c o / 



-0/(2-/3) 



f n hQ 
V Uh 



1/(2-/3) 



(28) 



In order to find the plasma parameters in the course of com- 
pression one should solve Eqs. (23), (26) and (28) for S(l), 
n c {l), Tih{l) from I = Iq to smaller I. The magnetic field dissi- 
pates completely when the sheet width becomes comparable 
with the box width, 5 ~ At the last stage, one is unable to 
separate unambiguously the total plasma volume onto a hot 
current sheet and a cold magnetized plasma therefore the 
presented equations may be only used to estimate roughly 
the final plasma parameters. Taking this into account, one 
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can simplify the problem even more and solve these equa- 
tions in the limit i5 < (, when they are formally applicable, 
and take the limit 8 = I in the obtained solutions. 

Let us introduce the dimensionless variable A = 8/1. 
In the zeroth order in small A, both Eq.(23) and Eq.(26) 
reduce to the same equation 



n c l = n C Qlo, 



(29) 



so the system is nearly degenerate. In order to find the sec- 
ond equation, one should expand Eqs.(23) and (26) to the 
first order in A and eliminate the zeroth order term. In- 
troducing one more dimensionless variable Y = n h /n c , one 
gets 



(30) 



Transforming Eq.(28) to the dimensionless variables A and 
Y and substituting n c from Eq.(29), one obtains 



A_ = /Y l \ 
A V Yl ) 



1/(2-/3) 



(31) 



Now one can eliminate I from Eqs.(30) and (31) to get finally 
AdY 



(Y- 



1 ¥^A +y + /3 - 4 = °- 



(32) 



The solution to this equation is 



/ A_\ 4 "' 3 = Y ( Y + S3 - 4 
VAcJ Y I Y + f3~4 



3-0 



At A > A (but still A < 1) Y goes to a constant, Y = 
4-/3, independently of the initial conditions. Taking into 
account that in the pulsar wind Y = 3 (LK), one can take 
for estimates Y = Yq = 3. 

Now one can estimate the compression factor, k = lo/l, 
necessary for the magnetic field to dissipate completely. Sub- 
stituting A = 1 into Eq.(31), one gets 

k = A%- 2 . (33) 

Substituting n c = kn c o into Eq.(27) yields an estimate for 
the maximal Lorentz factor attained when the magnetic field 
dissipates completely 



1 
Ao" 



(2-/3) 
208-1)' 



1/(2-/3) 



(34) 



where a = fe 2 n c o/(47rm e c 2 ) is the initial magnetization pa- 
rameter. It follows from Eqs.(16), (18) and (19) that in the 
wind upstream the shock A ~ 1/a therefore the particles 
may be accelerated up to large energies, j m ^> 1, even at a 
moderately large a. 



5 PARTICLE ACCELERATION AT THE 
TERMINATION SHOCK IN A STRIPED 
WIND 

At the termination shock, the flow is compressed and the en- 
ergy of the alternating magnetic field is released. Close to the 
equator of the flow, the average magnetic field is small and 
nearly all the Poynting flux is transferred to the particles. It 
was shown in sect. 2 that the downstream velocity is non- 
relativistic in this case, therefore the proper density of the 
plasma increases ~ F w times where F w is the wind Lorentz 



factor upstream the shock. For typical parameters (see sect. 
3) this compression factor significantly exceeds that given 
by Eq. (33); this confirms the conclusion that the alternat- 
ing magnetic field annihilates completely at the termination 
shock. 

The flow is decelerated in the shock by the pressure of 
the hot downstream plasma and magnetic field. In the col- 
lisionless shock the deceleration scale is about the Larmor 
radius of those particles, which make a major contribution 
to the downstream pressure. These particles penetrate up- 
stream by their Larmour radius and exert, via the magnetic 
field, the decelerating force on the upstream flow. In the 
standard MHD shock, the downstream temperature is about 
the particle kinetic energy in the upstream flow therefore the 
upstream particles penetrate about all the shock width im- 
mediately after they enter the shock. So there is only one 
characteristic spatial scale in this case, namely that of the 
Larmor radius corresponding to the upstream kinetic energy. 
In the shock in the striped wind, the particles gain energy 
from the annihilating magnetic field so that the downstream 
pressure is determined by particles with the energy signifi- 
cantly exceeding the upstream kinetic energy. The Larmor 
radius of these particles significantly exceeds not only the 
Larmor radius corresponding to the kinetic energy in the up- 
stream flow but even the strip width. Therefore the width of 
such a shock, or the deceleration scale, exceeds all "internal" 
scales in the upstream flow and hence the the field annihila- 
tion proceeds locally in the proper frame like in the plasma 
smoothly compressed by an external force. Only when the 
alternating field dissipates completely, the particle Larmor 
radius, calculated with account of both thermal and kinetic 
energy, becomes comparable with the shock width and the 
flow decelerates further on like in the standard shock. So one 
can roughly separate the shock into two zones. In the first 
zone, the flow is decelerated and compressed by the pressure 
of high energy particles entering from the second zone. The 
magnetic field dissipates in the first zone roughly according 
to the simple picture outlined in Sect. 4. The plasma heated 
by the field annihilation in the first zone enters the second 
one therefore the second zone resembles the standard shock 
with a hot upstream plasma. 

The compression factor, k, necessary for the magnetic 
field to dissipate completely was estimated in the previ- 
ous section. The continuity equation in the relativistic flow, 
ror = const, implies that the flow Lorentz factor at the 
end of the field dissipation stage is T d ~ F w /k. Adopting 
the outlined in the previous section picture of the particle 
acceleration by driven reconnection, one concludes that a 
power-law particle distribution is formed at this stage. In 
the proper plasma frame moving with the Lorentz factor 
r<j, this distribution extends from 7 ~ 1 up to 7 ~ j m . Now 
the alternating field is already dissipated and the flow de- 
celerates further on, from V ~ Td down to T ~ 1, like in 
the standard MHD shock. Applying the particle and energy 
flux conservation, nTv = const and eY 2 v = const, one finds 
the maximal energy in the particle distribution downstream 
the shock, -fmax ~ r<i7m. The upper limit on the minimal 
energy may be estimated from the condition that the en- 
ergy of the most of particles is simply randomized but does 
not changes considerably; then 7 m i„ ~ IV Making use of 



© 0000 RAS, MNRAS 000, 000-000 



8 Y.E.Lyubarsky 



Eqs.(33) and (34), one gets 

/ min ~ A ^Tu,; (35) 

7— ~ T^r ^^ry^J ' (36) 

Substituting /? = 1.5 and parameters of the Crab pulsar 
wind upstream the termination shock (see sect. 3), one gets 
7min ~ 600; ^max ~ 10 6 , which is roughly compatible with 
the parameters inferred from the observed spectrum. 

These simple estimates show that particle acceleration 
at the shock in a striped wind may form such a particle dis- 
tribution that the energy of most of the particles is signifi- 
cantly less than the upstream particle kinetic energy whereas 
the plasma energy density is dominated by a relatively small 
amount of high energy particles. These high energy parti- 
cles may be accelerated further on by the 1-st order Fermi 
mechanism thus forming a high energy tail at 7 > fmax- 
Recent investigations (Bednarz & Ostrowski (1998); Gallant 
& Achterberg (1999); Kirk et al. (2000); Achterberg et al. 
(2001)) have shown that in ultra-relativistic shocks, the tail 
is formed with the power-law index j3 — 2.2 — 2.3 compati- 
ble with the observed X-ray spectra of the Crab and other 
plerions. So the observed broken power-law spectrum with 
a flat low frequency part may be attributed to the parti- 
cle acceleration at the termination shock in a striped pulsar 
wind. 



The proposed mechanism may also explain why no low- 
frequency cutoff is observed in the Crab radio spectrum. 
At the kinetic energy dominated shocks, the power-law tail 
is formed only at the energies exceeding the downstream 
temperature, which is about the particle kinetic energy up- 
stream the shock. Therefore if the Crab pulsar wind were ki- 
netic energy dominated, only rather low Lorentz factor of the 
wind would be compatible with the radio data, which would 
imply extremely highly mass loaded wind. In the striped 
wind, most of the energy is contained in the magnetic field 
therefore the wind Lorentz factor is lower than in the ki- 
netic energy dominated wind. Moreover the presented here 
simple model of the shock in a striped wind predicts that 
in this case the downstream temperature may be consider- 
ably less than the upstream particle kinetic energy. There- 
fore the lack of the low-frequency turnover in the observed 
plerion spectra may be naturally explained within the scope 
of the proposed model. Of course the outlined here qualita- 
tive picture of the particle acceleration in plerions may be 
considered only as preliminary; it may be justified only by 
numerical simulations of the shock in a striped wind. 
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6 CONCLUSIONS 

The observed spectra of plerions from the radio to the 
gamma-ray band imply a very wide electron energy distri- 
bution, from less than few hundreds MeV to ~ 10 16 eV. 
Most of electrons are accumulated at the low-energy end of 
this distribution. Although the synchrotron life time of these 
electrons exceeds the plerion age, there is strong evidence to 
suggest that they are accelerated now together with high 
energy electrons responsible for the hard radiation from the 
nebula (Bietenholz & Kronberg 1992; Gallant & Tuffs 1999, 
2002; Bietenholtz et al. 2001). The observed flat radio spec- 
trum implies that the acceleration mechanism transfers most 
of the available energy to a small fraction of particles and 
retains most of particles at relatively low energy. 

It is proposed in this article that the flat energy distri- 
bution is formed in the course of the particle acceleration 
by driven reconnection of the alternating magnetic field at 
the pulsar wind termination shock. It is widely believed that 
just upstream the termination shock the magnetic energy is 
negligible as compared with the plasma kinetic energy be- 
cause plasma dynamics in the nebula implies small magne- 
tization downstream the shock (Rees & Gunn 1974; Kennel 
& Coroniti 1984; Emmering & Chevalier 1987; Begelman & 
Li 1992). However the plasma magnetization and dynamics 
depend only on the average upstream magnetic field so the 
upstream flow may be Poynting dominated provided most of 
the magnetic energy is associated with the alternating mag- 
netic field. In this case the average particle energy grows 
significantly at the shock where the magnetic field annihi- 
lates. Therefore the particle distribution with the power-law 
index (3 < 2 is formed readily. 
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